Implementing a topological quantum model with cavity lattice 
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Kitaev model has both Abelian and non-Abelian anyonic excitations. It can act as a starting 
point for topological quantum computation. However, this model Hamiltonian is difficult to im- 
plement in natural condensed matter systems. Here we propose a quantum simulation scheme by 
constructing the Kitaev model Hamiltonian in a lattice of coupled cavities with embedded A-type 
three-level atoms. In this scheme, several parameters are tunable, for example, via external laser 
fields. Importantly, our scheme is based on currently existing technologies and it provides a feasible 
way of realizing the Kitaev model to explore topological excitations. 

PACS numbers: 



I. INTRODUCTION 

Topological quantum systems are attracting consider- 
able interest because of their fundamental importance in 
diverse areas of physics, ranging from quantum-field the- 
ory to semiconductor physics 1 ^, as well as emerging fields 
such as quantum computation 3,4 . A topological quantum 
system has topological phases of matter that are insen- 
sitive to local perturbations. For example, a recently 
proposed quantum model with quasi-local spin interac- 
tions possesses emergent topologically ordered states 5 -. 
It has also been suggested that a two-dimensional (2D) 
electron gas in the fractional quantum Hall regime has 
non-Abelian anyonic excitations near particular magnetic 
fields, which can be used for topological quantum com- 
puting^. 

In addition to the search for naturally existing topolog- 
ical phases, there have also been wide-ranging theoretical 
efforts to design model Hamiltonians and artificial lattice 
structures that possess nontrivial topological properties. 
One of the prominent candidates for topological quan- 
tum computation is the Kitaev model on a honeycomb 
lattice 8 , as shown in Fig. |TJ This model can be described 
by a completely anisotropic Hamiltonian 
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where J x , J y , and J z denote the coupling strengths of the 
x-, y-, and z-type bonds, respectively. Using Majorana 
fermion method 8 or Jordan- Wigner transformationally, 
the Kitaev model can be analytically solved. It has two 
phases: a band insulator phase and a topological gapless 
phase 8 -. In the insulator phase, which is equivalent to a 
toric code model for the same braiding statistics**, the 
elementary excitations are Abelian anyons. However, the 
vortices in the gapless phase are not well-defined because 
the rule of the braiding statistics is unclear in this case. 
An applied external magnetic field breaks the time rever- 



FIG. 1: Schematic diagram of the Kitaev model on a honey- 
comb lattice. The building block in our scheme (denoted as 
a dashed box) is explicitly shown in Fig. Ufa). 



sal symmetry in the Hamiltonian and opens a gap in the 
gapless phase. Now the vortices obey a well-defined non- 
Abelian anyonic statistics^. With possible non-Abelian 
anyonic excitations, the Kitaev model could play an im- 
portant role in demonstrating anyonic statistics and in 
implementing topological quantum computing. 



Due to its strong anisotropicity, the Kitaev model is 
unlikely to exist in a natural homogeneous physical sys- 
tem. An alternative but promising way to realize this 
model is to use artificial structures 1 ^, such as Joseph- 
son junction arrays 13 , optical lattice d 14 ' 15 , and polar 
molecule o 16 ' 17 . However, with a Josephson junction ar- 
ray, the magnetic term of the model Hamiltonian cannot 
be cancelled completely^, while with other proposals (in 
either optical lattice or molecules) extremely strict con- 
ditions are required, such as ultra-low temperatures. 
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Recently, various theoretical proposals of using coupled 
cavities to simulate basic spin-interactions and many- 
body models have been studied^— . In particular, the 
XY spin model and the Hciscnbcrg model have been theo- 
retically simulated in reference^ and reference^, respec- 
tively. In reference^, the a x a x and a v a v couplings were 
proposed to implement by using two laser fields (here de- 
noted as A and B), with the related coupling strengths 
tuned by varying the applied laser fields. In addition, 
the a z a z coupling was proposed to implement by using 
another two laser fields (here denoted as C and D). How- 
ever, these three types of spin couplings cannot occur si- 
multaneously. To solve this problem, the Suzuki- Trotter 
formula was employed in reference^, which involves ap- 
plying an appropriate repeated sequence of pulses that 
tunes on the laser fields A and B for a short interval 
of time after another interval of time in which the laser 
fields C and D are tuned on. This can in principle yield 
an anisotropic Heisenberg Hamiltonian with these three 
types of spin couplings. 

In this paper, we propose an alternative approach to 
simultaneously implement all three different couplings of 
spins in a quantum model. However, in sharp contrast 
to the proposal in reference^, the simultaneous imple- 
mentation of these three different spin couplings in our 
considered model is simply owing to the introduction 
of suitable cavity modes depending on the geometrical 
structure of the cavity lattice. Moreover, our proposal 
is also experimentally accessible with currently existing 
technologies^—. Specifically, we show that a lattice 
of tunnel-coupled cavities can be tuned to emulate the 
strongly anisotropic Kitaev lattice. In our case, the basic 
element is a A-type three-level atom inside three cavities 
that are oriented in different directions, where the two 
long-lived atomic levels form an effective spin-i (qubit). 
The interactions between nearest-neighbor atoms are re- 
alized via the exchange of virtual photons between the 
coupled cavities. Indeed, our scheme for building qubits 
and achieving inter-qubit interactions is quite general, 
so that complex many-body models can be constructed 
by changing the structure of the cavity lattice and/or 
by varying the driving laser fields, without making use 
of the Suzuki- Trotter formula that involves applying an 
appropriate repeated sequence of pulses^. 

The paper is organized as follows. In Sec. 2, we de- 
scribe the design of our cavity lattice to simulate the 
Kitaev model. We derive the effective Hamiltonian of 
the system in two different cases in Sec. 3 and then dis- 
cuss the conditions for experimental implementation in 
Sec. 4. Finally, we give a summary of our scheme and 
our conclusions in Sec. 5. 



II. MODEL 

We construct an artificial 2D honeycomb lattice using 
tunnel-coupled cavities to implement the Kitaev model 
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FIG. 2: (a) Schematic diagram of the building block for con- 
structing the Kitaev lattice on a honeycomb lattice. At each 
site of the honeycomb lattice, there are three cavities oriented 
120° apart, where a A-type three-level atom is placed in the 
common region of the three cavity modes belonging to differ- 
ent cavities, (b) Schematic diagram of a A-type three-level 
atom coupled with two laser fields (denoted as gray arrows) 
and a mode of the cavity (denoted as black arrows) in the s- 
bond direction of the honeycomb lattice, where s = x and y. 
(c) Schematic diagram of a A-type three-level atom coupled 
with two laser fields and a mode of the cavity in the z-bond 
direction of the honeycomb lattice. In both (b) and (c), a 
laser field of frequency v\ drives the transitions |o) 44 \e) and 
\b) 44 |e) with Rabi frequencies Q, a \ and Qbi> respectively. 
The detuning between ui and ui ea = e e — e a ((jJeb = e e — £b) is 
A a i (A(,i). Another laser field with frequency v-2 also drives 
the transitions \a) 44 |e) and \b) 44 |e) with Rabi frequencies 
fi a 2 and Qb2, respectively. The detuning between 1/2 and 0J ea 
(u) e b) is A a 2 (A^)- Furthermore, a cavity mode of frequency 
vu is coupled to the transitions \a) 44 |e) and \b) 44 |e) with 
coupling strengths g% and g£, respectively, where k = x, y 
and z. The detuning between vu and uj ea {uj e b) is S% (Sb)- 



(see Figs. Q] and [2]). The building block of this lat- 
tice involves four lattice sites, as shown in the dashed 
box in Fig. [TJ At each site, we place an identical A- 
type three-level atom, which is coupled to three cavities 
[see Fig. Ufa)]. The three- level atom has two long-lived 
ground states, \a) and |6), which are denoted as spin- 
down and spin- up states, respectively. These effective 
spin-i's form qubits. The a x a x , a v a y , and a z cr z cou- 
plings are achieved by exchanging virtual photons be- 
tween nearest-neighbor cavities in the three bond direc- 
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tions of the Kitaev model. 

The 2D lattice of cavities with embedded atoms is de- 
scribed by a Hamiltonian consisting of three parts: Ha, 
He, and Hac- Here Ha is the Hamiltonian of the bare 
atoms: 



H A = 



N 

E 



(2) 



where a vq — \p)(q\, with p,q <E {a, 6, e}, j is the lattice 
index, and co ea {uJba) is the energy difference between 
states |e) (\b)) and \a). In Eq. ([2]), the energy of state \a) 
is chosen to be the zero point, and we set h = 1. 

He is the Hamiltonian of the cavity lattice in the ab- 
sence of atoms. In contrast to early worksi^r— , our model 
involves more cavity modes, and thus He can be de- 
scribed as 



N 

^ = E E 4 fc } (i)+ E Hi cl 

j k—x,y,z k—x,y,z 

r(fc) 



h cj = E ^(4*^ + °*4i) 



(3a) 

(3b) 
(3c) 



where ajy, with k = x, y or z, is the annihilation op- 
erator of photons, with frequency v k , for a mode of the 
cavity at the jth site, which is oriented in the fc-bond di- 
rection. For simplicity, we assume that all cavities in 
the same bond direction have the identical mode fre- 
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quency. H c ^(j) is the Hamiltonian for cavity photons 

at the jth site of the lattice, and H^j describes the tun- 
nelling between nearest-neighbor cavities along the k di- 
rection, with tunnelling strength tk- Here the small in- 
teractions among photons of different cavity modes are 
neglected, because the common region of the three cavity 
modes at each site can be designed very small. 

Lastly, Hac describes the interaction between atoms 
and photons: 



H AC = ^2[H^(j) + H^(j) 

3 

h^U)= E E (sNw4+h 

k—x, 

^ } (i)=E 
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where H^ t , with i — 1 (2), is the interaction Hamilto- 
nian between the atoms and the cavity photons (the two 
laser fields). Each cavity mode is coupled with transi- 
tions | a) «-» |e) and |6) O |e), and these transitions are 
also driven by two laser fields. In H^, is the cou- 
pling strength between the cavity mode k and the atomic 

(2) 

transition \a) <B- |e) (|6) o |e)). In H>J, Q al (f2 a2 ) and 



(^2) are the Rabi frequencies involving the driving 
processes of the transitions \a) o |e) and \b) <H> |e) by 
the laser field of frequency vi {1/2} ■ Here we assume that 

Uba > 9a, 9b, fiftl, ^a2, ^61, and ^62- 

In the following section, we show how to obtain the 
a x a x , a y <7 v , and <j z <t z couplings in the x, y, and z di- 
rections simultaneously, so as to implement the Kitaev 
model using the cavity lattice. 



III. EFFECTIVE HAMILTONIAN 

A. Reduced pseudo-spins 

For the atom at site j of the honeycomb lattice, when 
the interaction picture with respect to Ho(j) — Ha{j) + 

J2k=x,y,z H Cn(j) is USed > H inIU) and H intti) arC tranS " 

formed to 



E E( 
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where the detunings between the laser frequencies and 
the energy differences of the atomic transitions are given 
by A Q i = oj ea - vi, A a2 = uJea - v-i, Abi = uj eb - v\, 
and A{,2 = <^eb — ^2- The detunings between the cavity- 
mode frequencies and the energy differences of the atomic 
transitions are — u> ea — Vk and = u> e b — Vk- 
For a more detailed description, see Figures. 2(b) and 
2(c). In our scheme, all the detunings are red-shifted 
except for A&i. Here we assume that A i, A a 2, Aji, 
Afc2, <^a and $b are au i n the large-detuning regime, i.e., 
I^IJA^I > |^|,Oy3, where a — a and b, and j3 = ol, 
a2, 61, and 62. In this regime, two kinds of transitions 
can dominate when A i ~ 5^, A a 2 ~ 5 Z , Am ~ S^, and 
Afc2 ~ S^, where s = x and y. One is the Raman pro- 
cess, in which the atom is excited from state |a) (|6)) by 
absorbing a single laser photon and then falls to state 1 6) 
(| a)) by emitting a single cavity photon and vice versa. 
The other is the Rayleigh process, in which the atom is 
excited from state \a) (|6)) by absorbing a single laser 
photon and then returns to \a) (|6)) by emitting a single 
cavity photon and vice versa. Also, the Rayleigh process 
can occur by absorbing a single laser (cavity) photon and 
then emitting a single laser (cavity) photon. 

Because the detunings shown above are all in the large- 
detuning regime, the excited state can be adiabatically 
eliminated from the Hamiltonian as long as the initial 
state of the atom is within the subspace spanned by the 
two long-lived states \a) and |6). This yields the two 
long-lived states to be shifted in energy and also be ef- 
fectively coupled by the two-photon Raman and Rayleigh 
processes. We first study the simple case [ see Figs, ^h) 
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and 02c)]: 

A al =5 s b , A a2 = 5*, A b2 = 8% = (6) 

With the dominant transitions considered, Hj(j) + 
Hj 2 \j) is reduced, in the large detuning regime, to^ 

H e (j) = - (Val + Va2) Ola - (Vbl + Vb2)<J J bb 

- E (K"L"h<, ■ K"!/'^ 

k—x,y,z 



E [(M«+ 4 V aJ )a! j +H.c. 



- [(A zl al a + A z2 aQ a^+Hx. 
+ F(al l a n ) , 



where -qp = r2^/4A^, with (3 — al, <z2, 61, and 62; A* 
{g k a Y/5l and X b = (g^/8*; 

gffiai A _ gffi 
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A z2 = 
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g s b ^b2 



251 



Here we assume that X k = = X b . F(aj n a n ) in the last 
line of Eq. represents all the terms containing the 
operators (m, n £ {x,y,z}, and m 7^ n), and de- 

scribes the effective interaction between different cavity 
modes via atoms. Note that the term F(al n a n ) does not 
appear in the other proposals of cavity-based quantum 
simulations of spin models^—, because only one cavity 
mode was used therein. 

If 2\g a \ > then X k > Vfy/AAp. Thus, the AC 

Stark shift of two long-lived states induced by the laser 
fields can be neglected. As a matter of fact, when the 
condition 2\g a \ 3> is not satisfied, one can alter- 

natively introduce other laser fields to compensate the 
AC Stark effect^. In the large-detuning regime that we 
consider above, the zero-photon subspace is preserved be- 
cause only virtual two-photon processes are involved. As 
shown in the next section, all the photonic degrees of free- 
dom are eliminated in the zero-photon subspace when 
implementing the second adiabatic elimination. Thus, 
we can neglect the effective interaction between differ- 
ent cavity modes in Eq. ((?]), because after the second 
adiabatic elimination, the terms arising from F(al n a n ) 
can be guaranteed to vanish in the zero-photon subspace. 
To briefly summarize, in this subsection we have shown 
that the three-level atom at each site of the lattice can 
be reduced to an effective two-level system in the large- 
detuning regime. 



B. Effective couplings between pseudo-spins 

Below we attempt to eliminate the photonic de- 
grees of freedom and derive a pure spin Hamil- 
tonian for the pseudo-spins defined by the atomic 



states I a) and |6) of the original three-level atoms. 
We consider the interaction picture with respect to 

H = -J2k=x,y,z X k( a lj a kj^L ~ a lj a kj°bb)- When 

the tunnelling term between nearest-neighbor cavities, 
J2k=x,y,zJ2{i, j ) t k(al l a k:j +a ki a\-), is included, a new 
form of the Hamiltonian is obtained in this interaction 
picture. Again, using the assumption of the large detun- 
ing 
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we make the second adiabatic elimination^. Here the 
dominant process involves only the exchange of virtual 
photons between nearest-neighbor cavities (as shown in 
Fig. [3]), because other transitions, which oscillate rapidly 
with large frequencies, can be neglected in the long-time 
limit. Finally, the effective spin Hamiltonian in the zero- 
photon subspace can be written as 



H, 
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where s = x, and y. 



We define 7^ = 5 k /5 b , so 7^ = j y in this simple case. If 
the parameters are properly chosen as listed in Table. 1, 
the energy shifts of the two long-lived states, which gener- 
ate from the adiabatic elimination, can cancel with each 
other, and the coefficient of the magnetic term thereby 
becomes zero. Furthermore, these conditions lead to 
Jxi = Jx2, J y i = —Jy2, and J z \ = J z2 = -J z 3, so the 
effective Hamiltonian is reduced to the Kitaev model on 
a honeycomb lattice: 
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Cavity / Cavity j 



FIG. 3: Schematic diagram of the exchange of virtual photons 
between adjoining cavities to induce an effective coupling be- 
tween atoms at nearest-neighbor site of the honeycomb lat- 
tice. The gray dashed curve with arrows denotes the virtual- 
photon exchange. 



Table 1 Conditions of the parameters for imple- 
menting the Kitaev model in the case when Eq. ^ 
is satisfied. 
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are the x-,y- and z-type coupling strengths between 
atoms at nearest-neighbor sites when Eq. © is satis- 
fied. So far, the three different spin couplings in the 
Kitaev model have been emulated simultaneously by in- 
troducing three different cavity modes depending on the 
geometrical structure of the honeycomb lattice. 

Furthermore, by properly changing the Rabi frequency 
Sl a2 in the last column of Table. 1, a new magnetic term 
can be generated: 



C. More generic effective Hamiltonians 



In the above subsections, we have derived the effective 
Hamiltonians in the case shown in Figs. &h) and [He). 
In this subsection, we consider a more general case where 
the detunings for the driving laser fields are not equal to 
the related detunings for the cavity modes, i.e., A Q i ^ 5 b , 
A a2 ^ 6%, and A b2 ^ = 6%, as illustrated in Figure 4. 
The following new notations 
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where B = B z — (J zl — J Z 2)/4, and J zc = (J z \ + J Z 2 — 
2J Z 3)/4. With this effective Hamiltonian, it is possible 
to investigate the phase transition from a gapless phase 
of the Kitaev model to a gapped phase when the effective 
magnetic field is presented. 



5 x{y) = 5 x a M-A b2 = 5 x b iv) ~A aU S z ee S z -A a2 = S z b -A b2 

(12) 

will be used in the following derivations. We can show 
that the resulting effective Hamiltonian can be obtained 
via adiabatic eliminations without involving the unde- 
sired oscillating terms. In fact, under the new conditions 
in this case, all of these parameters Sk should be much 
smaller than the detunings for the laser fields and the 
cavity modes. 

Through two adiabatic elimination processes that are 
same as in the above subsections, we can obtain, within 
the zero-photon subspace, a similar effective Hamiltonian 
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e a , \a) 



e a , \a) 



FIG. 4: (a) Schematic diagram of a A-type three-level atom 
coupled with two laser fields and a mode of the cavity in 
the s-bond direction of the honeycomb lattice, where s — x 
and y. (b) Schematic diagram of a A-type three-level atom 
coupled with two laser fields and a mode of the cavity in the z- 
bond direction of the honeycomb lattice. In both (a) and (b), 
we introduce nonzero detuning differences S s = A a i — S§ — 
Aw — 5^ and S z = A a 2 — 5 z a = A b 2 — S z , in contrast to the 
case in Figures 2(b) and 2(c). 
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are the x-, y-, and z-type coupling strengths between 
atoms at nearest-neighbor sites when Eq. ([T^]) is satisfied. 

Similarly, an effective Hamiltonian containing the mag- 
netic term can be obtained in this more general case by 
changing the Rabi frequency fl a2 in the last column of 
Table 2, 

H eS = B'J2^ + J' x E 

j x— link 

+ J v E + E °l°h ( 15 ) 

y — link z — link 

where B' = B' z - (J' zl - J' z2 )/4, and J' zc = (J' zl + J' z2 - 
2^ 3 )/4- 



as in Eq. but with different parameters 
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where we have defined e = 1/(1 + Sk/Xk), k„ = 1 + 
4 = 1 + S k /2Sl 



S k /26%, and k\ = 1 + 5 k /25^. Now the condition (9) 



becomes 



a, + k >i^M( 1 + -L), tfc . ri.8, 



4 V^a A/3 



Table 2 Conditions of the parameters for imple- 
menting the Kitaev model in a more general case 
when Eq. (fT3)l is satisfied. 

Parameters x-type bond y-type bond 
{a x a x ) (a y ay) 



0& 



Parameters 



^62 
^a2 

06 



(g k a) 2 =ik(g k b ) 2 



z-type (<7 z a z 



{n z a n a2 f = lz (K z b n b2 ) 2 

9l9l < 



If the parameters are tuned to satisfy the prescribed 
conditions listed in Table. 2, we can also derive an effec- 
tive Hamiltonian for the Kitaev model, 

H eS = j' x E °i°* + Jy E °i°j+ J 'z E 
a; — link y— link z — link 

(14) 



So far, we have shown that the effective Hamiltonian 
of the cavity-lattice system with embedded atoms can be 
reduced to the Kitaev model if we choose proper detun- 
ings and Rabi frequencies for both laser fields and cavity 
modes. Note that in the cavity systems, every two ad- 
joining cavities can be linked by an optical fiber, and the 
photon-tunnelling rate t between nearest-neighbor cavi- 
ties can be continuously tuned by twisting the fibers. 
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IV. POSSIBLE EXPERIMENTAL 
IMPLEMENTATION 

In this section, we show how the scheme considered 
above can be achieved using the currently existing tech- 
nologies. Central to the success of our scheme is that the 
coefficients of the effective Hamiltonians must be much 
larger than the decay rates of the cavities and the excited 
states \ej). Below we show that these requirements are 
experimentally accessible. 

We use the following notations: ft — max{f2 al , f2 a2 , 
fi&i, ^2}, 9 = max{ ff *, g$}, A = min{|A al |, |A a2 |, 
|A(,i|, | A.52 1 } j S = rnin{<5fc}, and t = mhi{tk}. The occu- 
pation of the excited states \ej) can be estimated as [19] 

(l^fel) « (§) > (16) 
and the photon number n p h is 

The coupling constants between the pseudo-spins are ap- 
proximately given by 

(18 » 

The two types of the effective decay rates are represented 
by Tx = (n/A) 2 1 , and T 2 = \gSl/(SA + g 2 + n 2 /A)] 2 K , 
where 7 and k are the rate of the spontaneous emission 
from the excited level and that of the cavity decay, re- 
spectively. Hence, the condition ri,r 2 <C [<?Q/(£A + 
g 2 )] 2 t is experimentally realizable when 

7<minjt(J) j' and K<< - ( i9 ) 

Furthermore, the assumption (13) can also be satisfied 
in our scheme. By defining /1 = and r\ = |JA/pfi|, 

this inequality approximately reduces to fx + rj 3> 1/2, 
which can be satisfied if 2|<SA| > \gtt\ or 2\g\ > \Q\. 
In the ordinary case, detunings can be adjusted to be 
larger than gigahertz, while g and Q are around hun- 
dreds of megahertz for strongly coupled cavity-atom sys- 
tems. Then the first inequality, i.e., 2|<5A| ^> \gSl\, can 
thus be satisfied. On the other hand, for ultrastongly 
coupled cavity-atom systems, such as photonic band gap 
cavities^ 7 -, where g can reach ~ 20 GHz, the second in- 
equality, i.e., 2\g\ 3> can also be fulfilled. 

The above arguments suggest that cavities with a high 
cooperativity factor and a high cavity-atom coupling 
strength g are desirable for implementing the Kitaev 
model. Here the cooperativity describes the loss from 
the atom and the cavity, and is defined as C = g 2 /ivy. 
Encouragingly, cavities with the required properties can 



be realized in certain high-g microcavities with cur- 
rent micro-manufacturing and micro-etching technolo- 
gies (see reference^ and references therein). For exam- 
ple, in toroidal microcavities, the cooperativity factor of 
C ~ 10 7 was achieved experimentally 29 , and in the pho- 
tonic band gap cavities, C ~ 10 3 was also reahze d 25 i 27 ' 30 . 
In order to achieve the honeycomb lattice, one can build 
three optical cavities in different directions at each site of 
the lattice, or design an appropriate structure of the pho- 
tonic band gap cavities to restrict the photons of different 
frequencies in three directions. However, some microcav- 
ities have only one mode at each site of the cavity lattice, 
such as the toroidal cavities^, which are not good candi- 
dates for the implementation of the Kitaev model. Alter- 
natively, such cavity systems may be useful for simulating 
other spin models, e.g., Heisenberg spin models^—. In 
summary, we highlight that the high-finesse optical cav- 
ities (such as Fabry-Perot cavities )£L and photonic band 
gap cavitie o 24 ! 25 are two promising candidates to imple- 
ment our scheme in 2D lattice of microcavities, owing to 
the possibility of designing an appropriate geometrical 
structure for different cavity modes. 



V. DISCUSSION AND CONCLUSION 

Here we emphasize that the different couplings of the 
effective Hamiltonians can be simultaneously achieved by 
using suitable cavity modes depending on the geometrical 
structure of the considered cavity lattice and by properly 
applying detuned laser fields. Our scheme can also pro- 
vide a model system for observing the phase transition 
between the gapless phase and gapped phase of the Ki- 
taev model by tuning the external laser fields. In our 
proposal, for simplicity, all the parameters for the laser 
fields are assumed to be identical. However, if one can 
control the laser fields at each site of the honeycomb lat- 
tice, the model described in this paper can be realized 
with more flxibilities, owing to the tunable parameters 
of the system. 

In conclusion, we have proposed an approach to realize 
the Kitaev model by using an anisotropic cavity lattice 
with embedded field-controlled A-type three-level atoms. 
Our approach provides an artificial, but experimentally 
realizable, many-body spin system for demonstrating the 
topological phases in the Kitaev model. Also, our ap- 
proach can be extended to simulate other 2D and even 
3D spin models. 
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